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Abstract 

The issue of the observabihty of the Lamb shift in systems with non-degenerate energy 
levels is put to question. To this end, we compute the Lamb shift of such systems in the 
electromagnetic environment provided by two infinite parallel conducting plates, which 
is instrumental in demonstrating the existence of the so-called Casimir effect. A formula 
giving the relative change in the Lamb shift (as compared to the standard one in vacuum) 
is explicitly obtained for spherical semiconductor Quantum Dots (QD). It is the result 
of a careful mathematical treatment of divergences in the calculations involving distribu- 
tion theory, which also settles a controversy on two different expressions in the existing 
literature. It suggests a possibility of QD non-degenerate energy spectrum fine-tuning 
for experimental purposes as well as a Gedankenexperiment to observe the Lamb shift in 
spherical semiconductor quantum dots. 



PACS numbers : 12.20.Ds, 71.35.-y, 73.22.Dj 
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1 Introduction 



Both Lamb and Casimir effects were discovered or predicted in the late 1940s [HIl], and have 
been actively studied ever since [21 II]. These two phenomena are the most striking effects 
in standard Quantum Electrodynamics (QED) for they provide the strongest experimental 
support to the quantization of the electromagnetic field. 

In free space, the ground state of the quantized electromagnetic field is not just a empty 
state. It is the siege of quantum fluctuations which contribute to the so-called zero-point energy 
and lead to observable effects. In hydrogen-like atoms, the s-levels are the only ones getting 
dressed by the electromagnetic field quantum fluctuations and thereby are split from other 
states, which should have been degenerate with them according to Dirac theory, giving rise to 
the Lamb shift [H [6] . 

Quantum fluctuations also induce a Casimir effect. In general, the summation of the zero- 
point energy fluctuations of the electromagnetic field yields a divergent ground state energy. 
In the absence of coupling with gravity, this divergence is usually subtracted off in an additive 
renormalization scheme. However, a careful analysis on its volume dependence reveals the 
occurrence of a finite and observable force, known as Casimir force. The modification of the 
electromagnetic field boundary conditions in a fixed spatial volume, for example by the presence 
of two parallel perfectly conducting plates, should manifest itself as an attractive force between 
these plates [7]. Despite a previous (but not totally satisfactory) attempt at detecting the 
Casimir force [8], an experimental convincing proof of its existence has been only brought to 
light in the late 1990s for separation distances in the range of respectively 0.6-6/xm [9] and 
0.1-0.9/im [in]. In these experiments, one plate has been replaced by a perfectly conducting 
sphere. The replacement of one of the plates by a sphere allows to weaken the boundary effects 
due to the finite size of the plates, and thereby have facilitated its observation. 

In the context of Cavity Quantum Electrodynamics (CQED) [II1II2], our recent investiga- 
tion on cavity-induced effects on atomic radiative properties of quantum dots (QDs) has led us 
to wonder how much the Lamb shift of a system would change if put inside a Casimir device, an 
environment provided by two parallel conducting plates. The change of the atomic Lamb shift 
induced by a modification of the zero-point energy due to new electromagnetic field boundary 
conditions has been studied in [13]. It has been also shown how hermitian conditions are suffi- 
cient to separate the contributions of vacuum fluctuations from those of self-radiation reaction 
to the energy level shifts [H] . The coupling of an electric dipole to a conducting surface through 
absorption and emission of its own radiative field is a well-known problem from a classical as 
well as from a quantum- mechanical point of view [T^ [TB] . However, the understanding of the 
respective role of vacuum fluctuations and radiative field on the energy level shifts between two 
parallel plates seems to remain ambiguous [T7] . 

In this paper, to pursue investigations on the analogy between semiconductor QDs and real 
atoms, we propose to sketch a protocol based on the Lamb shift due to the modification of 
the zero-point energy to observe the atomic-like Lamb shift in spherical semiconductor QDs, 
established recently in [18]. This issue arises from the fact that, exception made of the azimuthal 
quantum number degeneracy for obvious symmetry reasons, the energy levels of an electron-hole 
pair confined in a spherical semiconductor QD are non-degenerate [19] . Then, the observability 
of Lamb shift in QDs is a conceptual problem, because the s- and p-levels are both shifted 
by the Lamb effect, as opposed to real atoms. This may explain the lack for theoretical 
works addressing the Lamb effect in such confinement structures. To this end, in section 
[2], expressing the change in Lamb shift for a system placed in a Casimir device is established 
by a careful mathematical treatment. In order to avoid masking the physics, the details on 
calculational steps are exposed in the appendices with the required rigor. The net outcome 
is beneficial because it settles a controversy between two different expressions found in the 
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literature [131 [T7]. Then in the following section [3l we apply the obtained results to the case of 
a spherical semiconductor QD, described in the formalism of the effective mass approximation 
(EMA) f2U\ fIT\ . The explicit dependence of additional Lamb shift on the separation distance 
of the plates in the Casimir device suggests a convenient way to modulate the QD energy 
spectrum when needed in experimental contexts. We then propose a Gedankenexperiment to 
discuss the possibility of observing the Lamb shift in semiconductor QDs, at least in a so-called 
strong confinement regime and for a judiciously chosen semiconductor. A concluding section 
summarizes our main results. 

2 Modification of the Lamb shift between Casimir plates 

Let us consider two parallel perfectly conducting squared plates of linear size L placed at a 
separation distance d ^ Lin vacuum. Intuitively, because the zero-point energy fluctuations are 
more important outside than inside the volume defined by the plates, a Casimir effect arises as 
an attractive force between them [2] . According to the phenomenological Welton argument for 
atoms [6] , the energy level Lamb shift may be viewed as due to the particle position fluctuations 
induced by the zero-point energy of the unrestricted surrounding electromagnetic field. Thus, 
when this surrounding electromagnetic environment is changed to that of a Casimir device, one 
should expect that the Lamb shift takes a different form and value for a given quantum state 
specified by a same set of quantum numbers. We argue that this difference in Lamb shifts may 
be exploited to make the Lamb shift observable in systems with non-degenerate energy-levels. 

The Lamb shift in real atoms placed in a Casimir device has been calculated with the Bethe 
approach for a relativistic electron in [13], as well as for a non-relativistic electronic motion in 
[17] . They both predict an additional shift to the standard Lamb shift, which depends on the 
separation distance d between the plates and goes to zero in the limit of d — )■ oo. One can 
reasonably expect a similar result in the Welton approach. 

However, there is a discrepancy in the leading contribution to the additional Lamb shift in 
these results [131 [T7|. For hydrogen- like atoms, it is shown in [13] that this additional shift is 
inversely proportional to the separation distance d, and that it is the sum of a non-relativitic 
contribution, scaling as the atomic Rydberg energy E*, and of a relativistic correction inversely 
proportional to the electronic bare mass rrie — in units with h = c = 1. In the non-relativistic 
limit, where the typical binding energy of the electron to the nucleus is negligible against its 
typical rest energy, i.e. if E* -C me, only the non-relativistic contribution survives, so that the 
additional shift finally scales as oc This is in flagrant contradiction with the result of [T7] . 
which predicts, in the non-relativistic limit, that the leading contribution scales as oc ^. In 
this work, we succeed to remove this discrepancy and show that the correct behavior goes as 
oc ^, bringing to light the reasoning mistake made in [T3] . 

2.1 Statement of the problem and assumptions 

Let us consider a theoretical non-relativistic spinless particle of mass m* and of charge ±ge, 
described by an Hamiltonian Hq, which is diagonal in some orthonormal basis of its eigenvec- 
tors {|n)}n with energy eigenvalues {En}n- The interaction of this particle with a quantized 
electromagnetic field is given by the Pauli-Fierz Hamiltonian in the Coulomb gauge [22]. In 
this section, we will make use of a convenient way to retrieve the standard Casimir force be- 
tween the plates, which is inspired by [23] but formulated and handled in the framework of 
distribution theory, in order to deduce the density of states modified by the Casimir device. 
The detailed and involved mathematical steps are presented in the two appendices [A] and [Bl 
This modification of the density of states induces the expected modification of the Lamb shift 
undergone by the energy levels of the considered particle. 
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To perform such calculations, it is only needed to evaluate the energy shift induced by 
the vacuum fluctuations, the effect of the particle self-radiation reflected by the plates being 
neglected. To this end, we choose a regime in which the separation distance d is sufficiently 
small to allow the emergence of the Casimir effect between the plates, i.e. d should be at most 
of the order of magnitude of /im [10] . Furthermore, we will assume that the typical wavelength, 
scaling as oc (i, due to the confinement between the Casimir plates should be greater than the 
order of magnitude of the wavelength associated with an authorized radiative transition between 
two energy levels of the particle. In such a weak coupling regime, the coupling of a two-level 
quantum atom with itself through absorption and emission of dipolar radiation reflected by the 
Casimir plates is dominated by the coupling of the two-level atom with the electromagnetic 
field vacuum fluctuations [151 [IS]- In a hydrogen-like atom, the associated Rydberg energy E* 
typically characterizes the radiative transitions, their wavelength scaling as oc Then, the 
weak coupling regime is fulfilled if we impose that 

Kd< K*, (1) 

where = ^ is the ground state energy in the presence of the Casimir plates. The quantity 
K* (X E*, being the usual Bethe average excitation energy used as the IR cut-off for the Lamb 
effect [SlEj, is surprisingly higher than the associated Rydberg energy E*, and though represents 
the maximal excitation energy the atom could access 

The coupling of the quantum particle placed between the Casimir plates with its own ra- 
diation field shall be also neglected when applied to semiconductor QDs. The validity of the 
condition Eq. ([1]) in spherical semiconductor QDs will be dealt in more details in section [31 
As we shall see in this particular case, the Bethe cut-off k* is also proportional but is of the 
same order of magnitude than the typical QD radiative transition energy. It is then also the 
natural candidate to depict the characteristic properties of the quantum system under study, as 
far as Lamb effect is concerned. Contrary to real atoms, wavefunctions of an electron-hole pair 
confined in a QD, described in the standard effective mass approximation (EMA), are restricted 
to the region of space defined by the QD boundary surface [191 [SQ]- Then the probability for 
the electron, and to a lesser extent for the hole, to tunnel from the inside part of the QD to its 
outside surrounding is vanishingly small. Actually, even if the confinement potential exerted 
on the electron-hole trapped in the semiconductor QD is more appropriately represented by a 
finite potential step [25] instead of an infinite potential well, this assumption remains reason- 
able, since the tunnel effect probability remains exponentially small. Thus, by construction, the 
probability of interaction between a photon emitted by the QD inside the Casimir plates and 
the QD itself is negligible in comparison to the probability of its interaction with an excitation 
of the surrounding quantized electromagnetic field. 

2.2 Standard Casimir effect 

As described in [23], the limit of perfectly conducting plates allows to consider the Casimir effect 
as a manifestation only of the electromagnetic field vacuum fluctuations by uncoupling non- 
ambiguously the zero-point energy of the electromagnetic field from the surrounding matter. 
The electric and magnetic fields E and B are now supposed to satisfy the continuity boundary 
conditions on the plates 

E|| = = 0, 

where the indices "||" and "_L" stand for the tangential and orthogonal components of the 
fields with respect to the plates. A convenient way to describe the Casimir effect consists in 
considering that the rectangular box B of volume V = L'^d, built from the Casimir plates, is a 
waveguide in the direction orthogonal to the plates, which shall be noted ^-direction from now 
on. Then, the set of TM and TE modes of the electromagnetic field is a natural functional 
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basi^. This means that, except for modes for which n± = 0, each mode of wave- number 
knun^ = A/kn + {Kdn±)'^ of quantum numbers {n»,n±) G x N* — the tangential wave vector 



is kii = ^nii and the orthogonal wave vector k± = Hdn± is a natural cut-off wave number 



due to the finite size of the waveguide in the ^-direction — has two possible polarizations. If 
periodic boundary conditions along the plates are also imposed, the zero-point energy of such 
electromagnetic field is then given by the divergent series 

E(L,d)= t„,„^+^^i|£= Y, 

(n||,nx)eZ2xN* nySZ^ (ny ,?ix)eZ2 xZ 

Reference |23] prescribes how it should be regularized. Recalling that even if the plates are 
supposed to be perfectly conducting, any conducting material is transparent to radiation at suf- 
ficiently high frequencies. Then, modes of arbitrary high frequencies do not actually contribute 
to the Casimir force between the plates. This can be described by introducing a dimensionless 

cut-off function k i — )■ where is a UV cut-off, in the expression of the zero-point 

energy E{L,d), as follows 

The function is assumed to verify 0(0) = 1, without loss of generality, so that, in the limit 
of a perfect conductor, i.e. k <^ k^, each term appearing in the sum defining the regularized 
Erf){L, d) goes to the term with the same quantum numbers of the sum defining non- regularized 
£'(L, d). It is very useful to suppose moreover that G iS(M), where iS(M) denotes the Schwartz 
space of rapidly decreasing functions on M. This assumption allows us to rigourously justify 
our calculation in the distribution sense, based on a generalized and corrected version of the 
one given in [13] . The details of the calculations are given in appendix [B], and lead to the 
following expression of the zero-point energy of the electromagnetic field in presence of the 
Casimir device 

d) = vUk f-J-) + v[ dk ^g(-) 4-) , (3) 

the standard Casimir energy between the two plates being 

7r2 L2 



Ecasirmr{L, d) 



720 d^ 



where the expression of the function g: s i — )■ arctantan7r(| — s)xr\z(s) is proved in appendix 

{1 if X G A 
, , . is the so-called characteristic function of 
otherwise 

the subset ACM. This calculation of the zero-point energy of the electromagnetic field shows 
how its density of states is modified by the presence of the plates. This modification can be 
viewed as a perturbation on the case of free-space, i.e. without Casimir plates. More precisely, 
the first term of Eq. ([3]) consists of the well-known contribution of the zero-point energy of 
the electromagnetic field in a region of free-space of volume V , because of its characteristic 
behavior as the third power of the mode eigenenergy k. The second term, which is the only 
one depending on the separation distance d, is then the contribution due to the presence of the 



^For more details, one can refer to sections 8.1 and 8.2 of [26]. In particular, there also exists a TEM mode, 
which is non-trivial only if the waveguide section is simply connected. 
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plates, denoted by -E^asimir(-^' particular, it is actually, as expected, independent from the 

regularization function 6 in the limit — — t- 0. This is coherent with the fact that the Casimir 
energy is a physical quantity, as opposed to the first term, which is explicitly regularized by it. 
Therefore, the second term of Eq. (|3]) suggests that the modification of the zero-point energy 
due to the Casimir device can be indeed entirely determined by a correction to the standard 



density of states p{k) = ^ in free-space [27], defined by 



i\k) = —g{ — 



(4) 



We insist on the fact that this function is not strictly a density, since it is not positive. However, 
this should be considered as a perturbation to the case of the absence of the plates in the sense 
that, for a fixed mode of energy k > and in the limit of infinite separation distance between 
the plates, the correction vanishes, and the free-space density of states p is recovered in this 
limit 



Pd{k) 



p{k) 



1^ 

TT k 



k 

9\ — 

f^d 



- 2k 



-> 0. 



This insures that the usual free-space properties should be retrieved in the limit d ^ oo. This 
is the reason why we will refer to the function pd as a density by language abuse. Furthermore, 
it is shown in appendix [B] that the computation is more general, and applies to any function F 
of the eigenmode energy k, when its mean value between the Casimir plates is to be evaluated. 
The total density of states between the two parallel plates is then the sum of the density of 
states p and a perturbation coming from the density of states pd, due to the plates. Therefore, 
to study the modification of a physical quantity due to the presence of the Casimir plates, we 
shall focus on the contribution due to the correction pd, given formally by the relation 



dkF{k)pdik)(f)\ 



the free-space contribution of the standard density of states p being set aside. The functions 
and as well as the indices (3 and 7 are defined in appendix [Bl intuitively they respectively 
contain the regular part of the function F under study, and the IR divergent part, to be 
appropriately regularized by appropriate physical arguments. This argumentation will allow to 
directly identify the modification of the Lamb shift coming from the presence of the plates in 
both Bethe [5] and Welton [6] approaches, as discussed in the two following subsections. 



2.3 Bethe et Welton approaches 

The Bethe approach to the Lamb effect is purely pertubative. The quantum second order time 
independent degenerate perturbation theory is applied to the Pauli-Fierz Hamiltonian Hpp, 
where the electromagnetic field is treated in the weak field limit [5]. Using renormalization 
arguments, in free-space, the Lamb shift undergone by any energy eigenstate |n) of the quantum 
system is found to be 

where m* is used as a natural UV cut-off, since the assumption of non-relativistic particle is 
assumed here. Historically, this predicts a Lamb shift for the hydrogen atom 2s-level, which is 
in excellent agreement with experimental values [5l |2l] . 

The Welton approach is also a perturbative approach, but is more phenomenological. It has 
the merit of giving a physical picture of the origin of the Lamb effect. More precisely, the Lamb 
shift is interpreted as a fluctuation effect on the particle position due to its interaction with 
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the surrounding electromagnetic field. These fluctuations Ar can be described as a continuous 
random variable, whose probability density is a three-dimensional centered isotropic Gaussian 
distribution of variance ((Ar)^) = ^^^log^, where k* is the Bethe IR cut-off, and m* is 
used as a natural UV cut-off consistent with non-relativistic assumption, discarding fluctuation 
modes of order of the particle Compton wavelength [6]. The particle then moves in a new 
effective potential (^(r + Ar)), averaged on the fluctuation distribution, whose flrst corrective 
term AV^(r) in the flne structure constant a is precisely the term giving rise to the Lamb shift 

AE„ = / d^r |(r|n)|^ Ay(r) = -^^^^ log ^(n|VV(r)|n). 
J 671 m K* 

2.4 Modification of the Lamb shift 

Eq. ([5]) is a regularized version of Eq. (5) of Bethe original paper |5] , written in the formalism 
we have introduced previously, and in the non-relativistic limit k* <^ m* 



™ /„,v72t..„m„n r' dfc p{k) 

k + K* k"^ 



Ai^n"'^ = ^^(n|VV(r)|n) / 



In the non-relativistic limit k* ^ m*, Eq. ([5]) is also a regularized version of Eq. (3) of Welton 
original paper [5] 

Ai?r°" = ^4 (n| VV(r)|n) %ik). 



3 m*^ Jf., k^ ' 

Therefore, invoking the remarks we have made in subsection \2.2\ in both cases, the correction 
to the Lamb shift due to the presence of the Casimir device should be evaluated by replacing the 
density of states p in absence of the device by the corrective term pd, and the divergent integral 
is regularized following the prescriptions of appendix [B] with the test function </> G iS(M). This 
reasoning is a shortcut to adapt Bethe or Welton original arguments to this new framework. 
Let us insist on the fact that, in free-space, both methods are equivalent and give the same 
results. In particular, they prescribe the same way to regularized the IR divergence of the 
previous integral. This will not be the same in presence inside a Casimir surrounding. However, 
considerations on the Lamb shift in the non-relativistic limit provide a natural UV cut-off when 
needed, being the mass m* of the particle under study, i.e. in this part we will set = m*. 

Now, considering Bethe approach, we recognize that f^*{s) oc fy{s) = ^, (3 = 1 and 

7 = — 1, where r)* = —. Whereas, considering Welton approach, we recognize that / oc 1, 
f-/{s) = ^, /3 = and 7 = —2. In both cases, theorem [U but not proposition [HI applies, and 
particularly, according to theorem 12. al — we have 7 G Z \ N — a second IR cut-off n^, = rj^, Hd 
is needed here. As far as Lamb effect is concerned, the Bethe IR cut-off n* is the perfect 
candidate, then we will set that = k*. Moreover, thanks to the weak fleld coupling k* > Kd, 
the regularization parameter 77* , deflned by k*, verifles that 77* = 77* > 1. Theorem 12. al then 
gives a Taylor expansion in powers of ^ of the modiflcation to the Lamb shift due to the 
Casimir plates for each prescription to regularize the wanted integral. 

By performing the Taylor expansion of the function s: s 1 — )■ f^{s)f^{sn*), this directly 
yields, for any r G N* 

^ES^id) _ i-imp+i) . J 1 



p=0 

where the function F is the standard Gamma function, and 

pWelton,n/ 



1 -1 ) I'^^nl ^1 
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These expansions do not converge in the meaning of power series, actually they are strongly 
divergent, and should be understood in the meaning of asymptotic serieso Asymptotic series 
are the typical kind of objets encountered in Quantum Electrodynamics. For example when 
applying the usual perturbation procedure or computing Feynmann diagrams at a fixed order 
the fine structure constant a, the obtained series makes sense at all fixed order of the 



m 



perturbation parameter a, but do not asymptotically converge as a power series of it. 
Let us remark thaP 



^-^Casimir i^) 



< 



AE. 



Welton,n 
Casimir 



id) 



AEn AEn 

This result may be expected, since the Bethe approach integrates renormalization arguments, 
on the contrary to Welton approach, which is more phenomenological. In Bethe approach, if 
we may say so, the logarithmic IR divergence is indeed regularized by construction, while in 
Welton approach it is regularized ad hoc. Then, when the usual density of states p is replaced by 
the correction due to the Casimir plates, this generates a new IR logarithmic divergence in the 
first case, and a divergence, scaling as oc ^ in the second. This explains why the modification 
is less important in Bethe approach than in Welton approach. 

While the expressions for the modification to the Lamb shift AE^'^^^{d) and AE^^^^^'^{d) 
differ, they show the same first order, which scales as oc The difference between the 
predictions made by Bethe and Welton approaches indeed becomes significant either if the 
two IR cut-offs and n* are close enough, or if we consider a sufficiently large order in 
so as to the expansion begins to diverge. We will therefore focus on value of rinW such 



that 



Iter I 



< 1. Under these assumption, orders higher than two in — are negligible in 



2r(2r-l)7?2'- 

comparison with the shared first order, which is then sufficient in the weak coupling limit 
K* > Hfi. Finally, when (n|V^y(r)|n) 7^ 0, i.e. when the energy level under study actually 
undergoes the Lamb effect, we deduce the relative modification to the Lamb shift due to the 
Casimir plates 



AE, 



Casimir 



{d) 



AEr, 



K* 



log-1 



> 0. 



(6) 



Because of the positive sign, this always actually gives rise to a enhancement of the Lamb shift 
undergone by the particle. This result calls for a physical explanation. When the separation 



^ While the expansions ^E^'^^^^^{d) or ^E'^^^^l^id) do not converge, it is possible to recognize the so-called 



^Welton, n / 
^Casimir ^ 

Stirling series, which is known to give the asymptotic behavior of the Gamma function in the neighborhood of 
\z\ — >■ 00. More precisely, from Eq. 6.1.42 [55], in the neighborhood of z — >■ 00, with | aigz\ < tt, for any r g N* 



E 



^ 2n(2n- 1)^2" 

n=l ^ ' 



logr(z)- (z- i)logz + z- ilog27r 
nr+i[z), 



where the asymptotic behavior of the remainder Rr{z) is given by 

\b2r\ 



Moreover, from Eq. 5.2.11.16 

for any n gW 



g(_i)Pr(p + i) 

p=0 



which implies in particular that 



\Rriz)\ < 



2r(2r- l)|z|2'-' 
Abel theorem for series [3D| and Stirling formula 6.1.39 



we deduce, 



1 

V2 



< 



r(n + ^) J_ 



n-l 

E- 

p=0 



-imp 



p!r(i) 



V2 



AE. 



Bcthc,n / 



,n/ j\ A T^Wclton.n/ 7\ 

Casimir ^ ^^Casimir 



AE„ 



AE„ 
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distance d decreases, the amplitude of the electromagnetic modes inside the Casimir plates 
increases, while their number is fixed. This leads to the reinforcement of the interaction of 
the quantum system with the quantized electromagnetic field, implying a strengthening of the 
Lamb effect. Moreover, this relative enhancement does not depend on the quantum state under 
consideration. However, it shows an explicit concurrence between the different scales of energies 
< <^ m*, and then of characteristic lengths of the problem, as we shall see in the next 
section. This re-summed expression is in agreement with [T7], because up to the smallest order 
in the dimensionless IR cut-off r]*, the correction AE^^^^^^^^^d) scales as oc ^ oc ^. Contrary to 
this reference, in the regime where the quantum system does not interact with its own radiative 
field, our approach allows to explicitly compute the proportionality factor. In particular, it has 
been possible to factorize the mean value of the Laplace operator of the potential (n|V^l^(r)|n), 
which is considered as a characteristic of the Lamb effect according to Welton approach. 

3 Observability of the Lamb shift in spherical semicon- 
ductor QDs 

An experimental protocol, according to which the energy levels dressed by the zero-point fluc- 
tuations energy with or without the Casimir plates are compared, should allow to overcome 
the need of degenerate energy levels, or of exactly computed energy levels. As known, the 
experimental observability of the Lamb shift in hydrogen-like atom is due to the s- and p-level 
degeneracy, when the principal quantum number is > 2, in the absence of interaction with 
the quantized electromagnetic fleld. The Lamb shift arises as a separation of the ns-spectral 
band from np-spectral band, while they should stay merged in absence of Lamb effect. So, how 
would an energy level Lamb shift be detected for quantum systems displaying no spectral band 
degeneracy such as a QD? In such systems each non-degenerate energy level is dressed by the 
quantum zero-point fluctuations of the electromagnetic fleld, forbidding the detection of the 
corresponding bare energy level. Eq. ([6]) gives the protocol of a Gedankenexperiment, which 
may allow to access the Lamb shift in semiconducting QDs, flrst calculated in |18j . 

The model we use to obtain this Lamb shift in QDs is an improved version of the standard 
EMA, where a pseudo-potential is introduced to partly remove the usual divergence of the QD 
ground state energy in such models with small radius. One electron and one hole, moving with 
their standard effective masses m* in the considered semiconductor, are conflned by an inflnite 
spherical potential well of radius R placed at the QD surface, and interact with each other 
through the Coulomb potential. The common approach to treat the interplay of the Coulomb 
interaction of the electron-hole pair, which scales as oc j^, and the quantum conflnement, which 
scales as oc is to use a variational procedure, for which two regimes of electron-hole pair 
should be singled out according to the ratio of the Bohr radius a* = of the exciton, fi being 
its reduced mass, to the QD radius R. First, in the strong conflnement regime, valid for a QD 
radius R < 2a*, the conflnement potential sufficiently affects the relative electron-hole motion, 
so that the interactive electron-hole pair states should then consists of uncorrelated electron and 
hole states. In the weak conflnement regime, valid for a QD radius R > 4a*, the electron-hole 
relative motion is quasi left unchanged by the conflnement potential, so that excitonic binding 
states appear, as if the electron-hole pair has not been conflned. However, the exciton has to 
be treated as a conflned quasi-particle of total mass M, and its center-of-mass motion should 
then be quantized. 

For this simple model, in the strong conflnement regime, it is proven that the predicted Lamb 
shift, in QDs of size experimentally synthetized and used, is of the same order of magnitude 
as the one in hydrogen-like atoms, at least for judiciously chosen semiconductors, such as for 
example InAs or GaAs, and then seems to be observable. Since this is the relative modiflcation 
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Figure 1: Modification of tlie Lamb sliift in splierical InAs microcrystals for d = l/im ( — ), 
0.5/im ( — ) or 0.25/im ( ). 




of the Lamb shift due to the presence of the Casimir device, the method of calculation of 
the Lamb shift in QDs [18] is not of interest. However, this provides the Bethe IR cut-offs 



c,h 



12m*. _R2 

c,h 



and K* 



^2^j^2 respectively for the electron, the hole and the exciton. 
Let introduce the electron and the hole reduced Compton wavelengths A* = :rr^ and the 



c.h 



radius R^-^ 



2 v s-^c.h- These are interpreted as the lower bound for the QD radius, allowing 
fluctuations of the charge carrier, according to the Welton approach, to be confined inside the 
QD. Then, it is supposed that R^^ < R < d, with the additional constraint < kI^{R) given 
by Eq. ([1]). In this context, we deduce from Eq. ([6]) the relative modification of the Lamb 
effect undergone by a semiconducting QD placed between the Casimir plates 



ae: 



ie,h 

Casimir 



AE, 



e,h 




R' 



■log' 



R 



valid for 



R^ 1 



(7) 



Then, there is a competition between the dimensionless ratios > 1 and 4 < 1, characterising 

c,h 

the problem under study in the strong confinement regime, which is characterized in turns by 
the ratio 4 < 2. Fi gure [1] shows the behavior of this modification inside InAs QDs for several 
values of the separation distance d. Only the electronic contribution to the Lamb shift is 
represented on this figure, because, in InAs, either the hole Lamb shift is negligible against the 
electronic contribution (heavy hole, ^ ^ ^ ^ 0.026) or both are almost equal (light hole, 

— ^ 0.41), the effective masses being themselves almost identical. 

mo ' ' 

Figure [T] indicates that if the separation distance d is chosen to be 0.5/im, the experimental 
observation of the Casimir effect is possible. The modification of the electron-hole Lamb shift 
between the Casimir plates is of about 5-10% in spherical InAs QDs of radius in the range of 
10-20nm, which is of reasonable experimental size [3T]. It is possible to enhance the amplitude 
of this modification. The first idea is reducing the separation distance d until the order of a 
few tenth parts //m, the problem being then that the radius R should be reduced accordingly 
to still satisfy the weak coupling regime condition. This effect alone almost leads to the total 
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modification of the Lamb shift in free-space of at least 50%, which seems significant enough to 
be observable. 

There finally exists other corrective effects, such as the reflectivity of the metal used, the 
roughness of the surfaces of plates and sphere, and the finite temperature, which has an impact 
on the Casimir effect, as described in [10]. By the same kind of reasoning as above, it seems 
also possible to account for them in our description. Since these are corrections scaling as oc ^ 
to the standard Casimir force, it is sufficient to consider the two first terms in Eqs. ([6]) or ([7]), 
which are the dominant terms of '^'^cagmi'-('^) g^^^^ j^g ^^.^^ correction, also scaling cis oc ^4 , is to 
be considered. 



4 Conclusion 

In this work, we have developed a comprehensive computation method to obtain mathematically 
rigorous results on the Lamb shift for non degenerate energy levels in semiconductor spherical 
quantum dots placed in a Casimir device. A significant deviation from the standard Lamb 
shift {i.e. in vacuum) found in [IB] is revealed. The explicit formula giving this deviation 
suggests a way to finely adjust the spectrum of a quantum dot for experimental purposes, 
by changing the distance between the Casimir plates, provided that such micro-mechanical 
operations are realizable. Moreover, the energy shift order magnitude, brought to light in this 
Gedankenexperiment, at least for judiciously chosen semiconductor and QD sizes, seems to be 
sufficient to indicate that such QDs may be appropriate systems to observe this modification, 
and that the Lamb shift in non-degenerate systems may not be out of reach. 

A On the evaluation of the sum of a trigonometric series 

In this appendix, we propose to compute the value of the series 

Esin(27rios) , , _ 
P 

It is well-known that these series are the Fourier series on M representing the function g: 
s I — )■ arctantan 7r(| — s)xr-^z{s). The function g is odd and 1-periodic on M, and is equal to 
g{s) = arctantan7r(| — s) = 7r(| — s) for s g]0, 1[. Since the function g is of class by parts on 
M, we deduce from Dirichlet theorem [3U] that the Fourier series of g is point-wise convergent 
on M, and that 

s(s) = J: i^^. V. e (8) 

pen* ^ 

The error made in [TB] is the statement that g{s) = tt{^ — s) for any s G M, from Eq. (A14) 
in appendix A of this reference, but, as seen abive, g is neither odd (because g{0) ^ 0) nor 
periodic. 

To be as most exhaustive as possible, let us prove rigourously the validity of Eq. ([8]), using 
the Euler-Maclaurin formula, as proposed in [13], which turns out to be really helpful to get 
a fundamental convergence result. This result will be established despite a cumbersome proof 
which involves tricky points of regularization theory of series and distribution theory. 



10 



A.l Notations and convergences 

For any fixed e > and any s G M, let the even functions gl G S(R) and be defined by 

9l{t) = ^ e--* and g^,{t) = Vt G M. 

Then, gl > g'^, point-wise on R. It is possible to justify directly this point-wise limit. Let 

e— 5>0+ 

s eR*, then 

/ (itglit) = ^signserf^|s| -> ^signs = / dt^°(t),| 

where the function sign = xr+ — XM- is the sign-function on M, and the function erf: x 1 — > 
-i= rf dte~*^ is the error-function on M, which verifies erf a; > 1 and lerfl < 1 on M. When 

s = 0, it is obvious that gQ = g^ = on R, and the previous expression is trivially satisfied. 
Consider now the two point-wise convergent series on R defined by 

r \ ^ sin(27rps) _ 2 n sin(27rps) 
9e{s) = 2^ g,{p) = 2^ e P and go{s) = g^^p) = . 



The first point of our reasoning is to prove the following proposition, given a simple conver- 
gence result of to go in the limit e — )■ 0"*". 

Proposition 1. 

g^ )■ 5^0; point-wise on R. 

Proof. For any e > 0, the function g^ is odd and 1-periodic on R, so that it is sufficient to 
restrict the study of the point- wise convergence on [0, 
For any e > and for s G {0, |}, we have g^is) = 0. 

Let e > 0. For any s g]0, |[, the series ^gl{p) are absolutely convergent — because 
|fl's(p)| < e"*^^, Vp G N*, where e"*^ g]0, 1[ — , and the series ^ ggip) are semi- convergent thanks 
to Abel theorem [30]. Denoting by ||-||^ = sup^ | ■ | the usual supremum norm on any subset 



C R, Abel theorem also yields, for any G N*, 



'Ep>N9siP) 



00 



< y 0. This 



implies that the series e 1 — > ^gl{p) are uniformly convergent on R+. Then, e 1 — > ge{s) G 
C°(R+), and in particular geis) > go(s). 

By parity and periodicity extensions, one obtains finally g^ > go, point- wise on R. | 

As we shall see, this point-wise convergence is not sufficient, and we would rather have 
a convergence of g^ to go in the limit e — )■ 0"'" in the sense of tempered distributions i5'(R), 
denoted as g^^ > go in 5'(R)0 The difficulty lies in the fact that the previous inequality 

e->0+ 

US'- (■5) 11^^ ^ \sinTTs\ ' ^^^^^ for ^riy s G R\Z and then almost surely on R, forbids the use of dom- 
inated convergence theorem, because s 1 — )■ j^it^] ^ '^ioc(^' '^^) = {h :M. — )■ R/ J^ds \h{s)\ < 
00, VK C R compact}. Therefore, it is not possible to directly deduce that g^ > go in 

P'(R). However, we prove in proposition |4] that there exists a constant C > such that 
< C, almost surely on [0, |], which solves the problem. 



^cf. Eqs. 2.5.3.12 and 2.5.36.6 

^For more details on theories of the Lebesgue integral and of distributions, one can respectively refer to [32] 



and [331 
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Proposition 2. 

g, y go m S'{R). 

Proof. By parity and periodicity and proposition ||5f.(s) < C, almost surely on M, then 
C0 e £^(M, ds) = {/i : M+ — >R/ J^^ ds \h{s)\ < cx)}, for any test function G S(R). Hence, 
the functions g^ are tempered distributions on M, for any e > 0, and the dominated convergence 
theorem applies, leading to the expected convergence 

/ ds(7,(s)0(s) — ^ / ds (7o(s)0(s), V0 G 5(M) ^ (7e ^ (7o in 



A. 2 Euler-Maclaurin formula and consequences 

The Euler-Maclaurin summation formula is an important tool in analysis. It provides an 
estimation of the sum 'Ylip=Q^ip) the integral J^^jdi h{t), h being a sufficiently regular 
function on [0, A^] with G N. Let us assume, for convenience, that h G C°°(M), then for any 
G N and any r G N*, we have 

f;%)= / dth{t) + ^^^^ + + E {^^^"-^Hiv) - h^'-'^m + R^'w, 

where the remainder is expressed as 

J[0,N] l^^J' 

Here, the function i?2n is the unique 1-periodic function which coincide on [0, 1] with the 
Bernoulli polynomial B2n, and the real numbers {62n}n6N* are the non-vanishing Bernoulli 
numbers, defined by = -B2n(0) = -62^(1)5 for any n E N* ^U\. The useful Fourier series 
representation stands (c/. Eq. 6.22 p. 1032 [M] ) 

B2r{t) = 2{-iy~\2r)\ y ^Q^^^f ^ VtGM and Vr G N*. 

nSN* ^ ^ 

Let s g]0, |] be fixed for the rest of the section until further notice, the case s = being 
trivial. Let e > be also fixed. 

Let us consider gl G 5(M). Since gl is even on M, ^^2n-i 9lit)\^^Q = 0, for any n G N*, and 
gl{0) = 27rs, the Euler-Maclaurin formula is, for any G N and any r G N* 

p=l ^[0,^] n=i ^ 

where the remainder is given by 

The mistake in [13] lies on the fact that both limits r, A^ — >■ 00 are taken without justifications 
for any s G M*. For s = 0, this limit yields a vanishing remainder. However, for s G M*, care 
must be exercised, since the dominated convergent theorem does not actually apply in this 
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case. As we shall see, it is easy to prove that it is possible to take the limit — t- oo. But, the 
limit Rr{s) of the remainder R^''^{s) does not go to zero in the limit r — oo, for any s G M, 
even after having taken the limit e — ?■ 0"*". 



Proposition 3. 



The limit N ^ oo of Eq. (i<'e,r,N) exists, does not depend on r E N* , and yields 

g,{s)= I dtgl{t)-'Ks + R\{s). (★,) 



Proof. Let r G N*. 

j2n-l 



First, gl G implies that gl{N), ^gl{t)\ V 0, Vn G N*. 

Second, gl G d/x), where the set Vt is either M+ or N*, fitted with its natural measure 
/X, being respectively the Lebesgue measure or the so-called Dirac comb A = Xlpez ^P' being 
the Dirac measure at p G Z. Then, by the dominated convergence theorem, /j^^jq d/if?^ > 

dfigl, so that 

J29liP)^-^Y.9liP) and / dtgl{t)—-^ [ dtgl{t) 
Third, \\B2r\\^ = \b2r\ < oo and gl G S'{R), hence B2r^gl G £i(M,dt), and 

Then, the limit iV — oo is well-defined, and yields, Vr G N* 

^7,(s)= / dtgl{t)-Tis + Rl{s)= [ dtgl{t)-7rs + Rl{s). 



Proposition 4. 

r/ie /zmzi e — > 0^ of Eq. (^l) exists, and 

g{s) = go{s), VsG[0,i]. (*) 
Moreover, there exists C > 0, such that \\g-{^)\\^^^^ < C, Vs G [0, |], and g^ > g in 5'(M). 



Proof. First, we recall that it is proved in proposition [T] that the limit e — ?■ 0+ of the quantities 
ge{s) and dtglit) are defined, their respective limits being g^^s) and |signs. Only the limit 
e — 7- 0"*" of the remainder R\{s) is left to be justified. Since the Fourier series representation of 
the function B2 are normally convergent on M, we can write 

After two integration by parts, for any e > and any n eN*, since s G [0, ^], we get 



TV \ TV TV 

KM = —< erf^fn + s) — erf — si 
2 \ 
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Then, for any n G N* and for any e > 0, we have the following (in)equalities 



(n+s) 




< hl{n) = V^r die"* < 

Jj'fn-s) V e / 1 

^ ,2 




e f 



1 

2 




7re 



(n-l) 



Then, we deduce that the function e i — > hl{n) is continuous on M+ by continuous extension 

2 

with h^{n) = — ^s(^) is continuous on Ml, and < hl{n) < vre"^^""^'' > 0. 



Second, we have ||/is(n)||E'"'^^ < vre ^\ Vn G N*, implying that the series e i — > ^ hl{ 

Finally, e 



are normally convergent on [0,1], and ||-R^(s) < Tr^^^j^e '^^'^ = 

Rl{s) G C0([0,1]), and Rl{s) > EneN- = 0- The limit e -^'o+ of Eq. is 

therefore well-defined, and yields 



9o[s) = vr( 



g{s), VsG]0,i]. 



Since (^o and g : x i — > arctantan7r(| — s)xM.^zix) are both 1-periodic odd functions on M and 
coincide on [0, |], they are equal on mH 

We have also proved here the existence of the constant C > needed by proposition [21 For 
any e > 0, J^dtgl{t) = f sign s erf and J^dtg°^{t) = f signs, then Wj^^dtg, 
So that, for any s G [0, ^ 



[0,1] 



2J' 



dtg-sit) 



[0,1] 



7r 



2tt 



1 -e- 



C. 



Finally, to sum up, we have just proved that g^, y g, point-wise in M, and in 5'(M). 



^In the case s €]0, i], we have proved that the remainder Rl{s) goes to zero in the hmit e 0+. This 
becomes not true as soon as we consider s € M*, but the same reasoning can be done to determine the new 
hmit. Using the same notations, we have that for any n £ N* , e > and s € M* 



Ktin) ~ — sign s< erf — =(n + Isl) — sign (?i — Isherf — = m, — Isl 
2 I Ve Ve 



> -sign s{l - sign (n - |s|)} = /i°(n). 



>o+ 2 



Now, for any n > \s\, we have that ^)||oo' ^ ^ ^ R,CcLSoning els cibovs, wg conclude thcit 6 

Rlis) e C°([0, 1]), for any s G M*, and 



Rlis) 



J2 /.°(n)==E/i^(n) = 7r 



if s ^ Z 
if s e 



sign s ;c „ ^ 7^ 0; 



where the function [•]: x i — [x] is the integer part function on M.. The limit e — > 0+ of Eq. (-k^) is stih 
weh-defined, and yields 



go(s) = TT 



^ + [s]-.s if s(^Z 







if s G Z 



Once again, the functions g and go a-re equal on M, because they are both odd and 1-periodic on K, and coincide 
with the function x i — > 7r(i — x) on [0, 
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B Mathematics of the Casimir effect 



B.l Density of states between Casimir plates 

In this appendix, we present in details the general result on the density of states in vacuum or 
between the Casimir plates. To this end, let F be a non-zero function. The quantity of interest 

is totally formal, where the two functions and satisfy the following assumptions. 

i. s I — y f{i^ds) is a dimensioned function of class C°° on of dimension /c"'^, where 
/3 > is defined by the asymptotic behavior f{k) = 0{k~^) for k — )■ oo. 

ii. s I — y s'^, with 7 G M, is a dimensionless function of class at least on R*. 

Let us remark that (/3,7) is uniquely determined. Moreover, the function G C°°(]R+) and its 
successive derivatives are bounded on R. 

B.2 Integral representation of the quantity H^{L,d) 

Euler-Maclaurin or Poisson formulas [30] , being the same equation once written in the tempered 
distributions formalism, cannot be used in the previous expression, because the function s± 1 — > 
F[KiiS^\±^ 4>(^^s\\±^, is not differentiable at s± = 0, when sy = 0, with = y/sy + s^. 
Let be < 5 < ?7, and consider the regularized expression 



For S|| E and s± G M, denoting for convenience s i = ^/s?, + + 5^, the regularized function 



(S||,sx) I — > ^^^^$^,(4^) belongs to S(M? x M) — the function s 1 — > f^{s)<p{^s) is of 

class C°° and rapidly decreasing on R \ [—77,77], for any 7/ > 0. So that, Euler-Maclaurin and 
Poisson formulas hold. 

In the limit L ^ d, this allows first to replace the sum over the quantum numbers ny G 7^ 
in the previous expression as an integral over with the dimensionless measure ^^^d^ky, 
where the equality has to be understood in the meaning of asymptotic series|§ The Poisson 



"For simplicity, consider the function G: k 1 — > ^ /,,,,._ fi ( V^'^ + i'^dV)'^ ) '/'( — ~ ) i ^-nd denote 

Gi: s I — !> ^-^j^^y{sri), with the notation s±s = Vs"^ ± S'^. Since the function G is even, the Euler-Maclaurin 
formula yields, for any r G N* 

Moreover, ^ J^dkG{k) = J^ds Gi{s), we then focus on the quantity 



2d 



(2r)! dt^"- \ L 



.2r \b2r\ ( d 
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formula is now applied to the function s± i — > — p — "^il-^l? )? for any sy G M , and yields, after 
making the change of variables ky — > ftrfSy 



$^ — A/Sit + n 



Lemma 1. 

The function (sn , n±) 



Sn + 71-5. + ^7^ ) integrable on (M^ x Z, d^sn 



dA(n_L)). r/ie sum J2nj_ez ^'^^ integral /jja d^sy m Eg. ('k^^s,!^^ can 5e inverted. 



Proof. For any n± G Z, we can writ^, after performing the change of variables s = ^^fsf+n^^+rf 



d^sii 



52 



<2vr||/^ 



< 27r||/^ 



< 27r||/^ 



< 27r||/^ 



ds s 



2 S- 



|<l>,(s) 



ds 



ds|<l>,(s)|. 



\%{s) 



where the dimensioned constant Cr = /jfds |G^'^''(s)| does not depend on L. This imphes that for any r € N*, 



the hmit 



can be taken, leading to 



riGZ 



L 



2d 



f dsGi(s) + 
Jr 



2r 



/jg dfc G{k) in the meaning of asymptotic series. 



Here, it is, in particular, not possible to show that 2'^'^ ( 4)^^ C2r > 0, for fixed L ^ d. In the 

neighborhood of r — >■ 00, using the asymptotic behavior of the Bernoulli numbers \b2r\ ^ 2 ('^^-jL , the asymptotic 

behavior of 2^''|^^ (■^)^' C2r is given by "^^^^^(2^)^^ ^2r- However, the successive derivative in the constant 
C2r will generally produce factors of the order at least of (2r)!, which dominate the geometric dependence in 

2^ in the remainder of Euler-Maclaurin formula for sufficiently large r. 
^7 



The last inequality is obtained by integration by parts 
ds 



\<S>^{s)\=2J Jnl + rf -7j / dtVi\%{t)\+2 / ds dtVi\<i>^{t)\ 



<2y^\ dtVt\<P^{t)\+2 dsVJ / dtVi\%it)\ 



\n±\ 



ds\^^{s)\ 
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Furthermore, |$7(s)| > 0, then 3iVj_ G N*, such that < ^, Vs > iVj_. Then, for 

any n_L G Z such that > Nj_, we have 



d^sii 



Sn + 



52 



52 



< 



2vr||/''| 



2 ' 



which shows that the series ^ Jj^a d^sy 
by Fubini-ToneUi theorem 



$7(^y/sJ + r22_ +772 



are convergent, and 



( JsJ + ni + i^) ^ ^ 



/ d^S|| (g) dA(n_L) 

f^Usl + nl + 6A 

+ n2 + 52 VV II 



i.e. (s||,n_L) 



< 00, 



$^(^^s2 + n2 +7^2^ G /:i(M2 X Z,d2s|| ®dA(nx)), where d2s|| 



dA(n_L) is the usual product mesure on R2 x Z. This leads to the inversion result by Fubini 
theorem. I 



Corollary. The function 6 — )■ H^^^^^^L, d) is continuous on [0, rf\. 



Proof. The function 5 J' / ^ 

is deduced directly from the proof of lemma [Tl which suggests that for any Sn G M2, 



^y^s|~Tri5^H-^j is continuous on [0,?]]. The result 



/s2 + n2 + 52 VV II y 



< 



11/^11 



00 



< 



ll/'^ll; 



s2 + iV2 (s2 + ni + r/2)2 ni(s2 + Nl) 



where (s||,n_L) 1 — > n^{s\+N\) '^XK.-.]-N^,N^[{ni_) G £^1*^x^*5 ds||®dA(n_L)), and (sy, ra^) 



(s2 + r/2)-^ 



+ ra^ + 772 



X[-^x,^x]K) e X Z*,ds|| ®dAK)). 



Let now e > 0, and consider the fully regularized expression 



H 



d,ri,€ 



{L,d) 



^7+1 



ds 



■ ji^ 



■<l>7(4^)e''"^'^e-^^'' 



^ ^ 4>,5,r],eJ 



As we will see later, the limits 6, e 
not possible to take the limit 77 — )■ 0"*" in Eq. ( iK'^^s,r],e 
$^ will not satisfy the assumptions of proposition |6 



0"*^ in Eq. (ir^^s,r],e^ well-defined, whereas it will be 



) and then Eq. ("^f 5„), when the function 
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Lemma 2. 

The function (s|| , s±, p) 



^^^^^$^(4^)e2^^P"^e-^p' IS mtegrable on {R^ x R x Z, d^sn 



ds± ® dA(p)). And, the sum J2n^& 
be inverted. 



an 



d the integrals Jjg,2 d^sy and f^ds± in Eq. iif^^s,rj,e^ 



can 



Proof. For any Sn G M^, s_l G M and p G Z, we have that 



< 11/ 



$7(4 



where, by Fubini-Tonelh theorem, and after performing the change of variable s — )■ s,,. 



/ d^sy ® dsi_ ® dA(p) |$^(4^)| e"^IPl < — ^ / 
Jm^xRxz 1 — e J^2 



d^sii (g) ds 



< 



1 -e- 



||x ^7V'='f? 



with Cr, = /jg dss|$-y(sr,)| < 00. Then, (s||,s_L,p) ^■ 
Z, d^sii ® ds_L ® dA(p)), and the inversion property is obtained by Fubini theorem. 



Lemma 3. 

The family (H^)^, where H^: s 1 — y Xlpez^^^'^^^®"'^^ ' /^'^ '^''^V ^ > 0, convergent in S'{R), 
and goes to the Dirac comb A = A. 

Proof. Let the functions /le: s 1 — > -^h{^) and h: s 1 — )■ -^/vre"^'^*^^, which satisfies dsh{s) = 
1, such that > 5q in i5'(M)|^ We apply the same reasoning as in the proofs of propositions 

[iMlto the series H^: s 1 — y Ylp£Z^'^^^^'^^~'^^ ■ Using the Euler-Maclaurin formula, we have, for 
any s G M, 

H,{s) = j dte^'^^'e-''^ + 2 XI / cos 2n7rt e^^'^'e"^*' = ^ h,{n - s) > 0. 



^This is well-known that the Dirac comb A is equal to its Fourier transform A(s) = 'Ylp^z^^^'^^^ ™^ 5'(M). 
Since h g and h^{s) = di e^"^'**e~'^*^ , inversion theorem holds and yields e^'^^^ = j^^At e^^'^^^h^{t). 

Moreover, for any (j) e 5(M), by the dominated convergence theorem, 

ds/ie(s)(/)(s) = / ds /i(s)(/)(Vis) ^0(0)= / As 5a{s)(j){s) ^ h, > 5o in 5'(R). 

From this result and since {s,t) I — !■ K{t)4>{s) e /:i(R2,dsdt) we get 



<ls(j){s)i / dse"" (/)(s) = / <ith^{t)(t){t) > / At5o{t)<p{t) ^ 5o = l in 5'(R). 

Using the translation property of the Fourier transform, we have 5p{s) = e^^P'*, for any p G Z, and by linearity 

of the Fourier transform in iS'(M), we obtain that A(s) — J^pez'^'^^^^'^ ■ Finally, using Poisson formula for 
(f) e 5(R), we deduce 

/"dsA(s)0(s) = V0(ri) V<^(p) = /"dsA(s)0(s) ^ A = A in 5'(M). 
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Moreover, the function satisfies H^{s) < ^ [\/7 ("^^ "I" e "I" -'-) "I" l-^O' ^ Then, for any 
(j) G 5(M), s I — > e £^(M,ds), I.e. H, G 5'(M)/And, by Fubini-Tonelh theorem. 



ds (g) dA{n) h,{n - s)\(f){s)\ = / ds H,{s)\(j){s)\ < oo, 
z Jr 

we obtain that {s,n) i — y h^{n — s)(f){s) G £^(M x Z, ds (8) dA(r?,)). Then, by Fubini theorem, 

dsH^{s)(f){s) = / ds(g)dA{n)h^{n-s)(f){s) = / ds(8)dA(n)/;,^(s)0(s-ra) = / </)(s-n) 
Jrxz Jrxz Jr 



Since s i — )■ Xlnez*^('^ — n) E S(R), the hmit e — ?■ 0^ of the right- hand- side of this expression 
exists and we have 

[ ds H,{s)(j){s) yy"(j){n)= I dsA(s)(/>(s) ^ H, in 

Jr ^^0+ ^^0+ 



Proposition 5. 



Hls^AL^d) - — — ^ H^Q Q{L,d) = jL,d). 



5,e-S>0 



Proof. By lemma H we have {s±,p) i — > ^^^^^$^(4^)e2i^P"^e-^?'' G C\R x Z, ds^ ® dA(p)), 
for any e > and sy G M^, then it is possible to inverse the sum J2pe:Z integral /j^dsj. 



Hls,,AL,d) = L 



..7+1 



47r2 



d'sii / ds 



Moreover, by lemma [3l we have, for almost any (sn, s±) G x 



<^\\f 



M4 



3 MOO 



where $: si — ^ sf^f (s^ + 7 + l) + |$^(s)| G ds). This upper-bound almost surely 
on X M is independent from 6 g]0,?7] and is also valid for 5 = 0. Then, (s||,s_l) i — >• 

^^^^H^{s±)^A4'^) G £^(]R2 X M,d^S|| O ds_L), so that Fubini theorem yields, for any e > 

and any 5 G [0, 77] 

7+1 <. ri3f ||-L\ 

d) = / d^sii ® ds^LS^H,^s^)%is^l^), 



and the function 6 1 — )■ H^^^^^^L, d) is continuous [0, r/]. Then, the limit 5 — j- 0+ exists 
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But, for any 5 E [0,r/], since G C°(M) is bounded on M, (s||,s±) i — > ^^^^H,{s±) G 
X M). Then, by lemma O ^^^H,{s±) > /^(4^)A(sx) in ^'(M^ x M). Therefore, 

e— >0+ 

since (s||,s_l) i — > $7(4"^) G 5(M^,M), H^g {L,d) > i/f^ o(L,rf), where by properties of 

the Dirac comb A, we have 

Hls,,AL^ d) = / d^sii ® d.^^^-l^A(.^)<|.,(.f ). 

This imphes in particular that H^q {L,d) )■ iffg Q(L,(i) = H^^ (L,d), the reciprocal 

being directly given by the corollary of lemma [H | 



Theorem 1. For any 5 G [0,//], 



7+2 



1 + 



9{s) 



lis 



Proof. By lemma [21 the sum ^^g^ and the integral /j^dsj, can be inverted in Eq. {'k^^s,rj,eh 
yielding 



sin 27r ps 
( 



-ep- 



,-7+1 
L2^ 
TT 



1 + 



9e{S} 
ITS 



By the same reasoning as in the proof of proposition O, since by proposition HI Qe > 9 in 

5'(]R), i.e. the limit e — )■ 0^ of the right-hand-side of the previous expression exists 



7+1 



dss' 



9is) 



ITS 



$7(5^). 



Proposition 6. 

Assuming that 7 > — 1, then 



H^iL,d) = V- 



.7+2 



ds s 



9is) 



TTS 



f{s)f,is)<P 



1±. 



V / dk[p{k)+p,{k)]F{k)(P 
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Proof. There exists S > such that l^-y+^'S)! < ^, Vs > 5. Let K C M_,_ a compact set, for 
any 77 G K and s > S 



TTS 



<ll/ 



00 



< 



3 11/^1 



2 s2 



Moreover, (r/, s) 1 — > $7+1(5,,) G C°(K x [0,5]), and G £^[0, -^l, ds). Then, there exists 
Ck,s ^ and 5 > 0, such that, for any 77 G K and any s G [0, S] 



1 + 



TTS 



< 11/^ 



<ll/^ 



lAK)! 



l'^7+l(^r,)| + 

s 

Ck,S H ^ {Al'S) + /7(5'maxK)} 



Then, r] 1 — )■ H^^.^ G C°(R+), and the hmit of Eq. {ir^^rj,r]) well-defined 



Let us identify the modification to the quantity H^^{L,d) due to the Casimir plates, in 
comparison with the free-space, as 



H 



_F,Casimir 



{L,d) = V 



K 



7+2 



y/P -t- {Kd7]y 



Theorem [T] and proposition [6] justify rigourously the heuristic reasoning made in subsection 
12. 2[ If 7 < —1, it becomes impossible to define the limit 77 — )■ O"*". Then, it becomes reasonable 
to assume that 77 is actually lower-bounded by 77,, > 0, defined by a second physical IR cut-off 

= ?7*Krf, competiting with the IR cut-off k^, due to the presence of the Casimir plates. For 
example, as mentioned in the body of this article, as far as the Lamb effect is concerned, this 
second IR cut-off is the Bethe IR cut-off k* > k^. 

Theorem 12. al give the general expressions for the quantity if J'^^'^™"^(L, d) in the limit of 
perfect conductor, i.e. when — — )■ 0, restricting our study to 7 G Z and to the following 

particular choice for the function , depending also on the IR cut-off k*. 



A 



with r]^, 



where A is a fixed dimensioned constant. While theorem I2.bl gives the general expressions for 



the quantity H 



Theorem 2. a. 



F,Casimir 



{L,d), under the same assumptions. 



Assuming that 7 G Z \ N, for any r > 1, in the limit ^ — )■ 0, 



H 



-FjCasimir / 



^7+2 r r , J2(n-1) / i 
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Proof. Since 7 < —1, let us consider Eq. 



This expression should be seen as the action 



of the distribution gxK+ ^ S'(R.) on the test function s 1 — > /|^(s)/^(s^J</)(^s^J G 5(M). How- 



E 



ever, we remark that the action on S{M.) of its derivative {gxR+)' = 
5'(M) is simpler. 

Denoting s 1 — > = /^°° dt/^^^(t)/^(t^j0(fjt^j G C°°(M), the problem is that 

^i*,ri* ^ '5(]R), nevertheless $5,,?,, G C°°(M) is rapidly deacreasing in the neighborhood of 00 and 
is bounded on M. Since the Schwartz class iS(M) is dense for the supremum norm in this 
set of functions, hereafter denoted i3+(]R), it is common to extend the action of a particular 
element of iS'(M) on iS(M) to i3_|_(]R). This can be done here, since the distribution {qxr+Y 
actually only involves positive arguments of El Then, the action of {gXR+Y ^ S'{W) on 
^i,,ri, £ -^+(11^) has a meaning, and we get 



H. 



FjCasimir 



{L,d) 



7+2 



Since ^^,,71* G C°°(M), the Euler-Maclaurin formula can be applied to the previous expression 

for any r > 1. Furthermore, we have that $^^'^(0) = S^fl(^)fiM\s=o + ^(f^)' 
for any n > 1. And, in this case, by performing a Taylor expansion of the function t 1 — > 
fli^)fiM e C°°(M+), we also remark that ^^^"^"^^^(O) = 2(„_i)+/3-y ) • Moreover, by two 
integrations by parts, the remainder satisfies 



(2r)! 



< 



l^2(r+l)| 

(2(r + l))! 

|&2(r+l)| 

(2(r + l))! 

I^2(r+l)l 



^2r+l 



+ 



dt 



2(r+l){t) ^(2(r+l))^^^^| 



dt 



dt2-+l 



flit)f,it,J(t> 

1 



(2(r + l))! 



O 



2r+/3--y 



2r+/3-7 



(2(r + l))! vr/; 
Let us study the behavior of the dimensioned constant Cr{C,*, ?7*), which is given by 

J2r+1 



dt 



dt2^'+l 



Since t 



G C'^ 



and 



for any t > 0, by direct cal- 



culation, we prove that jj^F+i f^^it) /^{t^^J = 0( ^2r+/-7+i ) in the neighborhood of 00. However 



^As the vector space may be though as the topological dual space of the vector space there exists 

a classical way, based on the underlaying weak-topology and density properties in Banach spaces, to extend the 
action of distributions. 

Let be * e B+{R). Since S(R) is a dense subset of for the norm 



3, 3(*jv)iv e 



a so-called 



regularization sequence, such that 4*^? 



> 5*. Then, by dominated convergence theorem, 



JR 



ds(<?XR+)'(s)*iv(s) = Y *^(^') 

pGN* 



*(0) 



- / dsl'(s). 



As this regularized limit is well defined, this allows to extend the action of the distribution 
any element of i3+(M), by setting 



€ S'i 



to 
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2r + /3-7 + l > 2(r + l) > 4, i.e. t ^ ^ 
theorem, since G 5(]R), in the hmit — — > 0, this imphes that 



'^^TTT f-yitrj,) E £^(]R+,dt). Fiom dominated conveigence 



< 



dt 
dt 
dt 
dt 
dt 



dt^^'+l 
d^" 



^2r 



—t 



-t 



dt 
d 



H V* 



dt 



2r+ 



2r+l 



^2r 



2r+/3-7 



dt 



i2r 



tfnt)f,^,{tr, 



dt 



l2r+l 



dt2' 



—tr, 



oil] 



where we use V* ^'i{v*t) = ^-li^)^ t > 0. 

On the first hand, if /3 — 7 = 1, z.e. if /3 = and 7 = — 1, the proof of the theorem ends 
here, since the integral J^^dt ^^Tr+i /-i(^i) is finite, independently from 77^,, for any r G N*, 

then Cr{^^,r]^) = 

On the other hand, assuming that /3 — 7 < —1, for any n G [0, 2r + 1] and any t > 0, we 
have 



Therefore, it is straigtforward to show that 



< 



r(/3 + 2r + l) ^ 

r(/3) 



2r+l 



2^- r(/3 + 2r + l 



dt 



Finally, reasoning as above, since we assume /3 — 7 > 1, the integral dt fi{t)^f^{ti 



IS 



finite, independently from 77*, for any n G [0, 2r + 1], then Cr{^*,ri^,) = 0^ ^2r+/3-7 /, 
the expected result, in the limit ^ — 0. 

Theorem 2.b. 

Assuming that 7 G N, for any r > [|] + 2, m the limit ^ 0, 



which yields 



Proof. Part of the proof of theorem 12. al can be adapted here, the difference being that we 
start with Eq. (*0 o,o) instead of Eq. ("^^0,0)- accordingly replace the function by 
the function $^,,o- Let us give the estimation of the remainder of the Euler-Maclaurin formula 
written for any r > [f] ^ 



(2r)! ' 



< 



\h2r\ 

(2r)! 



ds 



i2r-l 



ds 



2r 



-flis)f,{s)<l> 



(2r)! 
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The behavior of the dimensioned constant Cr{^*) is given by 



ds 



ds 

2(r-l)+^-7 



ds 



max T(/3 + n) . . 2(r-l)-7 2r-l 
^ ^2r-l "g[0.2r-l] ( — ) ^ 



n=0 



ds 



d" 



ds 



-/^(s)(/)(s) 



By the same reasoning as above, the integral J^^ ds \ ^f'yis)(j){s) \ is finite, independently from 

K^, for any n G |0, 2r — 1]. Since 2(r — 1) — 7 > 1, for any r > [|] + 2, the reminder now 
satisfies 



(2r)! 



which yields the expected result, in the limit ^ — > 0. I 
B.3 Casimir energy 

Here, we are interested in the computation of the standard Casimir effect using the formalism 
developed in appendix [Bl This consists in the determination of the contribution to the electro- 
magnetic energy gievn by Eq. (|2]) due to the Casimir plates, i.e. the term depending on the 
function g, which is the correction to the density of states in vacuum 



where we have directly identified F{k) = |, then deduce = |, f-yis) = s^, /9 = and 7 = 2, 
and then apply theorem 12. b[ From this result written for any r > 3, we obtain the correct 
Casimir energy in the presence of the plates 



e: 



Casimir 



{L,d) 



'720 rf3 



o( 



Casimir 



md). 
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